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Abstract. We generalize eigen-phase preserving (EPP) supersymmetric 
(SUSY) transformations to TV > 2 channel Schrodinger equation with equal 
thresholds. It is established that EPP SUSY transformations exist only in the 
case of even number of channels, TV = 2M. A single EPP SUSY transformation 
provides an M(M — 1) + 2 parametric deformation of the matrix Hamiltonian 
without affecting eigen-phase shifts of the scattering matrix. 



1. Introduction 

In this paper we study TV channel radial Schrodinger equation with equal thresholds. 
Such equation may describe scattering of particles with internal structure, for instance, 
spin [1] H ||. Supersymmetric (SUSY) transformations allow analytical studies 
of Schrodinger equation with a wide class of interaction potentials [H El [2]. 
In particular, the inverse scattering problem [8j [9] for two-channel Schrodinger 
equation with equal thresholds may be treated by combined usage of single channel 
SUSY transformations [TO] and eigen-phase preserving (EPP) supersymmetric (SUSY) 
transformations |11) . These transformations conserve the eigenvalues of the scattering 
matrix and modify its eigenvectors (coupling between channels), in a contrast to 
phase-equivalent SUSY transformations which do not modify scattering matrix at 

aii nana nana. 

In [3] the two-channel neutron-proton potential was reproduced by a chain of 
SUSY transformations, where the coupled channel inverse scattering problem were 
decomposed into the fitting of the channel phase shifts [TtH [T5] and the fitting of the 
mixing between channels. The fitting of the mixing between channels was provided 
by the EPP SUSY transformations. 

This paper extends the two-channel EPP SUSY transformations to higher number 
of channels. 

The paper is organized as follows. We first fix our notations and recall basics 
of SUSY transformations [T71 HH1 US] • Given the explicit form of a second-order 
SUSY transformation operator we study the physical sector of SUSY transformations 
between real and symmetric Hamiltonians. We analyze the most general form of a 
second-order SUSY transformation for the case of mutually conjugated factorization 
energies. Then we discuss applications of SUSY transformations to the scattering 
problems and calculate how /S-matrix transforms. 
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We start section [3] re-examining the conservation of its eigenvalues in the 
two-channel case. There is the following asymptotic condition for EPP SUSY 
transformations. The term with first-order derivative in the operator of EPP SUSY 
transformation vanishes at large distances. 

To generalize EPP SUSY for arbitrary number of channels we study the matrix 
equation which comes from this asymptotic condition. We find that EPP SUSY 
transformations may exist for 2M-channels only and obtain their general form. A 4- 
channel example explicitly shows how EPP SUSY transformations acts. We conclude 
with a summary of the obtained results and discussions of possible applications. 



2. Second order SUSY transformations 

2.1. Definition of SUSY transformations 

SUSY transformations of stationary matrix Schrodinger equation are well known 
[TTl [19] . In this subsection we just fix our notations. 

Consider a family of matrix Hamiltonians EI = {H a } 

H a = -I N ~ + V a (r), (1) 

where Jjv is the N x N identity matrix, V a (r) is the N x N real symmetric matrix 
potential. A multi- index a parameterizes a family of potentials V a {r). A matrix 
Hamiltonian defines the system of ordinary differential equations 

H a ip a (k,r) = k 2 cp a (k,r). (2) 

on N x N matrix functions (p a {k, r). 

A (polinomial) SUSY transformation of Q is a map of solutions 

Lba ■ <P a {k, r) -> ip b (k, r) = L ba ip a (k, r) , (3) 

provided by a differential matrix operator 

d n d n ~ x d 

L ba = A n —+A n „ 1 — + ... + A 1 - + A , (4) 

where Aj, j = 0, ...,n, are some matrix valued functions. This differential matrix 
operator obeys the intertwinning relation 

LbaHa = HbLba ■ (5) 

The intertwinning relation ^ defines both the operator L ab and the transformed 
Hamiltonian Hi,. We will consider a family of Hamiltonians H[iJ a ] = {H \L oa H a = 
HbLta} related with the given Hamiltonian H a by transformation operator Q. In the 
next subsection we present explicit form of the second order transformation operator 
and the transformed Hamiltonian. 



2.2. Second-order SUSY algebra of matrix Schrodinger equation 

Given initial matrix Hamiltonian Ho, we choose two N x TV matrix solutions 

H Uj = EjUj , j = 1,2, (6) 

with Ei , Ei called factorization constants. These functions determine a second order 
operator L20 [1*1,1x2]: 

WW = [l N d 2 r -V0 + E1 + {E 2 - E x ){wi - w 2 )- 1 {wi - d r )} f(r) , (7) 
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where 



w j (r)=u' j (r)uj 1 ( r ), w\ + w' + Ej = V , .7 = 1,2. (8) 



Functions Wj(r) are called superpotentials. We also introduce a second-order 
superpotential 



W 2 (r) = (E 2 -E 1 )[w 1 {r)-w 2 {r)\ \ (9) 
The more symmetric and compact form of formula ([7| reads 



W(r) 



-H H h VK 2 d r 



f(r) ■ (10) 



Operator L 2 q and Hamiltonian i?o obey the following algebra 

L 2 qH = H 2 L 20 , H 2 = H — 2W2 , (11) 
Ll Q L 20 = (H Q - E X ){H - E 2 ) , L 20 L\ = (H 2 - E X ){H 2 - E 2 ) .(12) 

The new (transformed) potential is expressed in terms of second-order 
superpotential W 2 as follows 

V 2 = V - 2W2 . (13) 

The transformation operator has a global symmetry 

L 20 [ui,u 2 ] = L 2Q [ui(r)Ui,U2U 2 ] , detC/i i2 ^0. (14) 

The second order SUSY transformations of an initial Hamiltonian Hq form the 
family IH^-ffo] = {^2^20-^0 = H 2 L 2 q}. We will work only with Hamiltonians from 
HM-ffo] and we omit subscripts in the notation of transformation operator L 2 q — > L. 



2.3. Restrictions to the SUSY transformations 

We restrict our consideration only to second order SUSY transformations with 
mutually conjugated factorization energies E\ = E% = £. Potentials Vo and V 2 are 
supposed to be real and symmetric. Hence the transformation functions u\ and u 2 
have to be mutually conjugated, u\ — = u. The symmetry of V 2 demands the 
symmetry of superpotentials wj = w\ = w, w\ — w 2 = w* . Defining the 
Wronskian of two matrix functions as 

W[ui,u 2 ](r) = u\ \r)u' 2 (r) - u\' \r)u 2 (r) (15) 
= uj(r)[w 2 (r)-wf(r)]u 2 (r). (16) 

we see that the symmetry of superpotential w implies a vanishing self- Wronskian 
W[u,u] = of transformation functions [2]. 

We present the second-order superpotential W 2 in terms of the matrix Wronskian 
for further needs, 

W 2 (r) = [E x -^ 2 )u 2 (r)W[ Ul , M2 ]- 1 (r) U f(r). (17) 

To specify acceptable choice of transformation solutions explicitly, we choose the 
basis in the solution space. Natural basis for the radial problem, r £ (0, 00), is formed 
by the Jost solutions /(±fc,r) with the exponential asymptotic behavior 

/(fc,r^oo)^/jve ifcr . (18) 

Let us expand the transformation functions in the Jost basis 

u(r) = f (-K,r)C s + f (K,r)D , (19) 
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where K = k r + iki, K 2 = £, ki > 0. Complex constant matrices C and D should 
provide vanishing self-wronskian W[ii, u] — 0. The wronskian of two solutions with the 
same k is a constant. For instance, W[/(— k, r), f(k, r)] = 2ikl^- Then, calculating 
W[u,u] we get a constraint on the possible choice of matrices C and D, 

D T C = C T D . (20) 



Matrices C and D have an ambiguity due to symmetry (14). Rank of matrix 
C, rankC = M < N , determines the structure of transformation operator. The sum 
of ranks rankC + rankZ? > N, otherwise operator L is undefined. Using (14) we 
may transform C to the form, where only first M columns are non-zero and linearly 
independent. Reordering channels (by permutations of rows in the system of equations 
([2])) we can put nontrivial M x M minor of C into the upper left corner. Then, C and 
D obey the following canonical form, 

c =U'o)- »=(U°1)- 

where X = X T is a symmetric MxM complex matrix, and Q is (N—M) x M complex 
matrix. This canonical form is a gauge which fixes ambiguity ( 14 ) of transformation 
solutions. 



2.4. Application to the scattering theory 

In concrete physical applications of SUSY transformations we may further restrict 
the class of Hamiltonians. In particular, in scattering theory [T] we work with the 
radial problem, r £ (0, 00). The interaction potentials decrease sufficiently fast at 
large distances and may contain centrifugal term 

lim r 2 V(r) = 1(1 + I N ) , I = diag(Z x , . . . , l N ) , e ll7T = ±I N . (22) 

The physical solution has the following asymptotic behavior 

ip(k, r -> 00) oc /c~ 1/2 [e- ifer e y 5 - e ikr e~ il ^ S(k)] , (23) 

where matrix coefficient S(k) is the scattering matrix. 
Scattering matrix is related with the Jost matrix 

S(k) = e a ^F{-k)F- 1 (k)e a ^ , (24) 

where the Jost matrix reads 

F(k) = lim [/ T (fc, ry] [(21/ - I)!!]" 1 . (25) 

r— fO 

Diagonal matrix v indicates the strength of the singularity in the potential near the 
origin 

U(r^0) = v{v + I N )r- 2 + 0{l). (26) 

Knowledge of the Jost solutions allows one to define scattering matrix. Super- 
symmetric transformations of Hamiltonian and solutions induce the transformation of 
scattering matrix. Formal approach to the calculations of the S-matrices was devel- 
oped in the work of Amado [5U] . 

Let us consider how the Jost solution transforms asymptotically, 

(L/o)(fc,r->oo)= (27) 



2 t E2 + El t I TXT W l + W 2 



-kf + — - + I W 2 ) (r -> 00) - ikW 2 {r -> 00) 



exp(ikr) . 
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Assume that there exists the following limit 

U oa (k) = lim \-k 2 + E2 + El 

Then the transformed Jost solution reads 

/ 2 (fc,r) = (L/ )(fc,r)f/ co 1 (fc), 



Wl + W 2 



ikW 2 



(28) 



(29) 



Making similar manipulations with the physical solution ( 23 1 we establish the 
form of transformed S-matrix 



S 2 (k) 



(30) 



In the case of our second-order SUSY transformation, the transformed S-matrix 
depends on the factorization energy £ and parameters Q, X through the matrix 
multipliers Uoo(k) and U 00 (k) . In general, this dependence may be very complicated. 
Moreover, the scattering matrix S 2 may have unphysical low and high energy behavior. 

SUSY transformations that deform the scattering matrix in a simple way are 
useful tools to solve inverse scattering problem. In the two-cannel case there is a 
special kind of deformation, when Uoo{k) becomes an orthogonal matrix [TT]. We call 
such deformations as eigen-phase preserving transformations. 



3. Eigen-phase preserving SUSY transformations 

3.1. Two channel case 

Let us analyze conditions that make a two-channel SUSY transformation be an eigen- 
phase preserving one In this case parameters of the transformation, Q = q, X = 
X, are just some numbers. Matrix U^k) depends on q only and becomes orthogonal 
when q = ±i. The determinant of u vanishes at large distances, detu(r — > oo) — > 
with such choice of q. Let detu(r — > oo) ~ e, then superpotential w diverges as 
wlr — > oo) ~ e _1 and two-fold superpotential W 2 vanishes as w 2 (r ~ > oo) ~ e. As a 
result, the limit (28) contains only even powers of k 

2 , E 2 



Uoo{k) = lim 



Wl + w 2 



(31) 



2 2 

The cancelation of odd powers of A; is a necessary condition to provide EPP SUSY 
transformations. In the next subsection we establish the most general form of matrix 
Q which leads to the vanishing limit 

lim W 2 = , (32) 

r— > oo 

for the case N > 2. 

Parameter x is also should be fixed to provide dctW[u, u*] ^ for all r > which 
leads to a finite V 2 . 



3.2. Asymptotic SUSY transformation at large distances for arbitrary N 



The transformation function u ( 19 ) has the following asymptotic behaviour at large 
distances 



u(r oo) Uoo {In + Ar 1 + o(r 1 )) , 



Ae 



iKrY. 



(33) 
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Im 

Q 



-Q T 
In-m 



J M,N—M 



hi 
—In-m 



(34) 



For each concrete EPP transformation N and N — M are fixed, therefore we will use 
notation E instead of Y<m,n-m- 

The two-fold superpotential behaves asymptotically as 

,* 1-X..T 



lim W 2 = W 2 ^ = (£- £*)<>W[t 



2i£ lm A*e tK * rl: W[u 



-iKrT, A T 



Using asymptotic ( 33 1 we see that this limit is a constant matrix 

-l 



W 2iOQ = 2iS lm A*e iK * r * [t& «,)' - (u T Ju 
= 2i£ lni A*e tK " rT -\e- lKrY -A T A* ' 
= 2£i m A* [K* A T A*Y> + KEA T A*] _1 A 



iK*rT.\i 



z- iK ^\'A T A*e iK * 



(35) 
(36) 

(37) 



We introduce auxiliary matrix W c 



Woo := K*A T A*T, + KT,A T A* = 2 



k r (I M + Q T Q* 
iki(Q-Q*) 



) iHQ T - (Q*) T ) 
-k r (QQ* + I N -m) 



Limit ( 32 ) leads to the following matrix equation 
^2,00 = =}► A*W^A T = , 



detWoo ^ , 



(38) 



(39) 



which provides asymptotic cancelation of fc in (28). In the two channel case these 
equations fix Q uniquely. When N > 2, these equations determine a set of Q values. 

Equation (391 may be satisfied if and only if matrix A is singular. Matrix W x 
is invertible, rankWoo = N. Let rank A — n, then rank A* = rankA T = n and 
rank(W^ c 1 ^4 T ) = dimImg(W^ 1 ^4 T ) = n. The dimension of kernels dimKerA = 
dimKer A* = dimKeryl T = N - n. Equation (39) implies that Img(W^ 1 A T ) C 
Kery4*, hence n < N — n. Therefore equation (39) has solutions only if n < hN. From 
the other hand, from explicit form of matrix A, ( 34 ) , its rank n > max(Af, N — M) . 
That is, (39) has solutions if and only if 

N 

rankA=y, N = 2M . (40) 



From here it follows that for odd number of channels equation (39) has no solutions. 
Consider 2M x 2M matrix A 

it 



A = 



Im -<-. 
Q Im 

with rank A — M. Two its rectangular sub matrices have the same rank 



rank 



Im 
Q 



rank 



Im 



M . 



(41) 



(42) 



We can take first M columns of A as linearly independent, then from (40), (41) and 

(43) 



(42) follows that there exists M x M matrix Z, such that 



Im 
Q 



IM 



M . 



Solving this equation we obtain Z = —Q T and QQ 1 



'M- 
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Q = ±iB , B T B = BB T = I M , (44) 



and substitute Q in this form into ( 39 1 . First of all we invert matrix Woo ■ This matrix 
can be factorized in two ways 

k r B T (B + B*) -ki(B T +Bi) \ _ f K(B^ +B T )B* -K(B T + B^) 
~K(B + B*) -k r B(B^+B T ) )~\ -ki(B + B*) -k r (B + S*)St 

k r B T -kiI M \ ( (B + B*) 
-kiI M -KB )\ (St + B T ) 

(St +B T ) \f k r B* -KI M 

{B + B*) J { -KI M -ArBt 



(46) 



We note that 

W 00 W (X5 -4(fc r +/e l ) ^ Q (_B + B*) J \ (fit + B T ) (45j 

Therefore the inverse matrix reads 

w -i = _L_ ( KB^ -kj M \((b t + st)- 1 o 

00 2|#| 2 ^ -KI M -KB* ) \ (B + B*)- 1 

1 / k r B^{B T + St)- 1 -fc 4 (S + B*)- 1 
2\K\ 2 \ -^(£ T + 5t)-i -k r B*{B + B*)- 1 

Let us introduce notations for auxiliary matrices 

(S T + St) o 

(B + B*) 



B = 



Bi 



k r B' —KIm 
-KI M -KB* 



Then W 2 ,oo = £i m A*S fc S-M T /(fc r 2 + kf). 
Using the following matrix identities 

B T (B T + B^y 1 = (B* + S) _1 S* , B\B T + B^y 1 = (B* + B)- X B , (47) 
(S T + S t )- 1 S T = S*(S*+S)- 1 , (S T + S t )- 1 S t = B(B* +B)- 1 , (48) 

B A \ (B + B*)- 1 ) { -W Im ) (49) 

I M iB* \((B T + St)- 1 \ 

-iSt I M )\ (B + B*)- 1 ) ' 

we see that 

A*B k B- 1 A T = 

K < , st -»/«• \ / / M ? ;s* \ / (s T + st)-i o )=(| 



iIm -B* )\-iB^ Im )\ (B + B*)- 1 
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Thus ( 44 1 gives solutions of ( 39 1 . 

Now we can calculate the asymptotic form of a transformation operator explicitly. 
To calculate asymptotic of W2W we use the symmetry of superpotential w = u'u^ 1 = 

(u T r\u T y, 



lim W2W = W2,, 
r— > 00 



)Wco = 2i£i m u* 00 W{iir 
-2iK£ lm A*W-}ZA T 



37 M Oo] U Oo( U CX>) ( U CX>) 

-if In 



2£ In 



-iBt 



hi 
-iB* 



K* 



-A*Bi.B~ HA = 



(B- 





-B* 



Matrix W2.00 is real, therefore 



2 /lo 



^(B T - fit) 
-2Jm 



= Re(W 2 ,, 

21m 

i(B — B* 



(B T + Bt)-! 





(B + B*)" 



(51) 

4k 2 k 2 I N , and 



The matrix C/oo defined in ( 28 1 reads 

C/oo(fc 2 ) = (-k 2 + k 2 - k 2 ) I N 

Matrix f2 is real, orthogonal (up to a normalization), £! T f2 
antisymmetric ft — —Q T . To establish its orthogonality and antisymmetry one should 
use relations (47 1, (48). With these two properties of fl the matrix U OD (k 2 ) becomes 
proportional to the orthogonal matrix 

U^ik^UKk 2 ) = ((-fc 2 + k 2 - k 2 ) 2 + 4k 2 kf) I N (52) 

That is the Jost solutions at large distances are rotated by orthogonal matrix Uoo 

(Lf)(k, r — > 00) — > Uoo exp(ikr) , (53) 



In this case the S-matrix transformation (301 is just an energy-dependent 
orthogonal transformation, 

R s (k 2 )So(k)R^(k 2 ), 



S 2 (k) 

with the orthogonal matrix, R^Rs = I 



Rs 



S — *Ni 

m u k 2 



[(- 



■„2\2 



4fc 2 fc?r 1/2 



(54) 



(55) 



That is we obtain desired generalization of two-channel EPP SUSY transformations. 

The above analysis is valid for an arbitrary M x M symmetric matrix X. 
Transformed .S-matrix 5 2 depends on matrix Q only. Therefore X might provide 
additional M(M + l)/2 parametric deformation of potential V2 without affecting the 
S-matrix. From the other hand, possibility of such deformations contradicts to the 
uniqueness of the inversion of the complete set of scattering data. Therefore, there may 
exist only one matrix X corresponding to one physical potential V%. The EPP SUSY 
transformation should be uniquely determined by the factorization energy, M x M 
complex orthogonal matrix B and a sign factor. In the next subsection we show how 
to fix matrix X and prove that the corresponding potential V2 is regular for all r > 0. 
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3.3. Eigen-phase preserving SUSY transformation near the origin 

To analyze the properties of EPP SUSY transformation in the vicinity of r — we 
will use the solution 

Mk,r) = ~ [f (-k,r)F (k) - f (k,r)F (-k)] , (56) 

vanishing at the origin 

Mk, r^O)^ diag [^ TWV - ' (2^ + 1)!! ) ' (57) 



where Fo(k) is the Jost matrix (25). We rewrite transformation solution in the basis 
(tp (K,r), f (K,r)) expressing f (-K,r) from ((56 1 

2K 

fo(-K,r) = — ^{K,r)F^\K) + h(K,r)F a {-K)F^{K), (58) 



and substituting in (19) 



u(r) - 2J ^MK,r)F - 1 (K)C + f (K, r)(D + s C) . (59) 



where 



so - F (-K)F -\K) = ( 81 52 , e y f s e"f = S . (60) 

We can transform matrix (Z) + soC) to the form of matrix D (without affecting C) 
multiplying u from the right 

(D + s C)( , M = f * (61) 

where 

1 = X + S! ±i(s^B + S T s 2 ) -B T s 3 B, (62) 
Then the transformation solution reads 

u(r) = 2 ^MK,r)F -\K)(lf B J ) +/o(#,r) ( f T *f ? ) . (63) 

Consider the case X = 0. In this case the potential V2 is regular for all 
r > 0. Let us prove this. According to the Wronskian representation of the 
second-order superpotential W 2 (16), the potential V 2 will be regular if and only 
if detW[u,u*](r) ^ 0. 

The derivative of the Wronskian W[u, u*] reads 

W[u,u*]'(r) = {£ -£*)u T (r)u*(r). (64) 

By construction W[it, it*] is an anti-Hermitian matrix, i.e. W[it, it*] — — W^[it, it*]. We 
represent transformation solution in a block-diagonal form 

u(r)=(y U12 V (65) 
y u 2 i u 22 J 

with M x M matrix blocks. When X = these blocks obey the following boundary 
conditions: itn(0) = 0, it2i(0) = 0, ui 2 (oo) = 0, u 22 (oo) — 0. As a result the 
Wronskian 

w M - ( 3 § ) ( ) ~ ( 3 3 ) ( 5 1 :l 2 ) , (66) 



a 



12 
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has vanishing blocks at the origin, Wn(0) = 0, and at infinity W^oo) = 0. Boundary 
behavior of W\2 is not determined. 

Now we can calculate diagonal blocks of the Wronskian integrating its derivative 



W[u,u*]'(r) 



W[ x W[ 2 
W'\ 2 Win 



> u 1\ u 1\ u \\ u \2 ' u 2\ u 22 
X * , T * T*iT* 
U 12 U 11 ' U 22 U 2\ U Y2 U \2 ' U T1 U T2 



Integration of (67) with the established boundary conditions yields 



W[«,u*](r) = {£-£*) 



V 



/ (wfiWu + W2W21) dt W\2 


00 

W\ 2 - J (ul 2 u* 12 + ul 2 u* 2 



dt 



(67) 



(68) 



/ 



Assume that there is a point ro where detW[u, u*](ro) — 0. Hence, matrix W[u, u*](ro) 
has at least one zero eigenvalue 

W[u,u*](r )v = 0. (69) 

Let us represent TV dimensional eigen-vector v as two M dimensional vectors v u and 
Vd and rewrite ( 69 ) as a system of equations 

W 12 v d = , 



W n v v 



(70) 

W{ 2 V u + W 22 ti d = Q. (71) 

The first term of the scalar product (v u ,WuV u ) + {v Ul W\2Vd) — 0, ((a, b) = a*l>>) is 
positive 

{v u ,W n v u ) = 

ro r 

(v u , / dt (liu^! + u 21 u 21 ) v u ) = I ((u^Vu, u* n v u ) + (u 21 v u ,u 21 v u )) dt > , 



therefore (v u , Wvz^d) — n u < is real and negative. Now calculating scalar product 
(v d , W\ 2 v u ) + (Hd, W 22 v u ) = with negative second term 

(Vd, W 2 2V d ) = 



<0, 



we obtain a contradiction, (t^, Wj 2 u u ) 



n u > 0. This contradiction proves that 



Wronskian W[it, w*] have only non-zero eigenvalues for all r > 0. As a result W[it, u*] 
is invertible, and hence both W2 and V2 are regular (finite) for all r > 0. Any non-zero 
A will lead to the potential V2 which is singular in some point ro. 

This prove completes our construction of multi-channel EPP SUSY transforma- 
tions. In the next subsection we present an illustrative example. 
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■ 1234*1234tl2!4 




012340123401234 



Figure 1. Off-diagonal entries of t he exactly solvable potential matrix V2 
obtained from the uncoupled potential | |72[ l with parameters ai = 1.1, 0,2 = 1.5, 
as = 2.1, (14 = 2.5, b r = 2.5, bi = 1.3, for three choices of the factorization energy 
£ = —2 + 1.5i; — 1.25 + 3.i;4.5i. The strength of coupling increases with arg£ 
decreasing from 0.78-7T to 7r/2. 




1 2 3 4 12 3 4 




Figure 2. Diagonal entries of the exactly solvable potential matrix V2 obtained 
from the uncoupled potential l |72| with parameters a\ = 1.1, 02 = 1.5, 03 = 2.1, 
04 = 2.5, b r = 2.5, bi = 1.3, for three choices of the factorization energy 
£ = -2 + 1.5i; -1.25 + 3.i; 4.5i. 



3.4- ^-channel coupled potential 

We construct our initial 4-channel potential with I = 0, 

V (r) = diag [v (r, ai),un(r, a 2 ),v (r, a 3 ),v (r, a 4 )] , (72) 
from four copies of the following single channel potential 

vo (r,a) = ; 2 ° . (73) 
smn (ar) 

Its scattering matrix is diagonal and reads 

S (k) = diag [s (k, ax), s (k, a 2 ), s (k, a 3 ),s (k, a 4 )] , s a (k, a) = ^— . (74) 
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Consider an ingoing wave in jth channel 

ipinj = exp(-ikr)(6ij , S 3j , 8 3j , &y) T , (75) 

where is the Kroneker delta-symbol. Ingoing wave is just a first term of long-rage 
asymptotic 

[ e -i*r _ e ikr S Q (k)] (8^,8^,8^,8^ . (76) 
Scattering of such wave on the potential ( 72 1 results just in a phase shift of the outgoing 
wave 

=exp(ikr + 2iSo(k ) a j ))(6 lj ,6 3 j ) 83 j ,84j) T , (77) 

by eigen-phase 

k 

8o(k, a,j) = — arctan — , (78) 
a j 

without mixing between channels. 

Using EPP SUSY transformation we deform potential to introduce coupling 
between channels. Diagonal components of the basis (<po, /o) explicitly reads 

<Po(k, r; a) = — ^ =■ (k cos(fcr) — a coth(ar) sin(Ax)) , (79) 

f n \ ft .k + iacothiar) 

fo(k, r; a) = exp(zfcr) 1 —. . 80) 

k + la 

These solutions together with matrix B depending on a single complex number 
b = b r + ih, 



fi — Sin 



(82) 



completely define EPP SUSY transformation. 

Let us fix all parameters of the model aj,b r ,bi except the factorization energy 
S (ai = 1.1, a 2 = 1.5, a 3 = 2.1, a 4 = 2.5, b r = 2.5, bi = 1.3). In Figures [l] 
and [2] we show the potential V2 provided by EPP SUSY transformations for three 
values of S = —2 + 1.5i; — 1.25 + 3.i;4.5i. The strength of coupling increases with 
arg£ decreasing from 0.787T to ir/2. One can check that arg£ = corresponds to 
zero-coupling, V2 = Vq. For our choice of matrix B and parameters, the matrix f2 
reads 

/ -0.936848 0.305791 -0.16973 \ 

0.936848 0.16973 0.305791 

-0.305791 -0.16973 0.936848 

\ 0.16973 -0.305791 -0.936848 J 

The matrix S2 ( 54 1 has the same eigenvalues, but non-diagonal character of 
potential results in the mixing of different channels in the outgoing wave, 

[ e - ifc V't - e lkr e- { ^R s Sv{k)R T s ] (6 lj ,8 aj ,6 3j ,5 4j ) T , (83) 
There is another set of ingoing waves 

=exp(-ifcr)i?,(fc 2 ), R s = {R lt R2,Ra,R4), (84) 
given by columns Rj(k 2 ) of matrix Rs which scatter just with a phase shift (78), 

VWj = exp(ifcr + 2\5 {k, aj))Rj(k 2 ) . (85) 

Vectors Rj depends on the energy of ingoing wave. In Figure [3] we show this 
dependence for a particular example S = 4.5i. Changing two complex parameters 
b and E we can manipulate transitions between channels which may open a way for 
broad physical application of EPP SUSY transformations. 
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4. Conclusion 

In the present paper we have generalized two-channel eigen-phase preserving SUSY 
transformations to the multichannel case, N = 2M > 2. It was surprising, that 
such generalization exists for even number of channels only. A single EPP SUSY 
transformation depends on a complex factorization energy £ , and M x M complex 
matrix B, such that B T B = Therefore single EPP SUSY transformation provide 
an M (M — 1) + 2 parametric deformation of scattering matrix without affecting eigen- 
phase shifts. 

There are several possible applications of presented results. One can use EPP 
SUSY transformations to solve inverse scattering problem by deforming a diagonal 
S-matrix as in [3J. We also may consider the S-matrix eigenvalues which conserved 
under M(M— 1)+2 parametric deformation as integrals of motions for some dynamical 
system associated with matrix Schrodinger equation [4 . In this context it is interesting 
to establish how this dynamical system looks. We expect that in this way new exactly 
solvable non-linear equations may be discovered. 
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